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ABSTRACT 



ll has been shown {VIedwin and Beaky, J. Acoust. Soc. Am., v. 86. 1 124-1130 
(1989)) that spilling breakers in the laboratory produce a sound close to the Knudsen 
sea surface noise spectrum from 400 IIz to 20 kHz. The surface spectral production 
density of newly-created near surface bubbles under these laboratory spilling breakers 
has now been acoustically determined by using an array of hydrophones. The surface 
spectral density has been obtained by identifying the individual bubbles that create this 
noise spectrum. The radii were calculated from the resonance frequencies. The bubble 
positions on the surface of the water were determined from the difference in time of ar- 
rival of the bubble radiation to the elements of two vertical arrays of hydrophones. The 
production area and rate of production of bubbles of radii 0.048 to 7.40 millimeters have 
been calculated and the total volume of air encapsulated into bubble foam per unit area 
of spilling breakers has been determined. 
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I. INTRODUCTION 



The purpose of this tliesis is to develop a method for the determination of the bubble 
production density from breaking waves. The motivation for this study came from the 
observation that the spectrum produced from the noise of the breaking waves in the 
laboratory tank at the Naval Postgraduate School closely resembled the Knudsen sea 
surface spectrum of the ambient noise at sea. This development suggested that the 
breaking waves experimentally produced will exhibit the same bubble characteristics as 
breaking waves at sea. It is hoped that the capability of comparing the calculated bub- 
ble production density with the area of foam coverage monitored with satellite imaging 
and aerial photography will aid in acoustical and meteorological research of air sea 
interactions by predicting the number and size of bubbles which may be found in the 
region being studied. 

The method used in this research is to utilize passive acoustics to identify individual 
bubbles created from breaking waves. The breaking of the wave, characterized by the 
initial formation of a large amplitude, low frequency bubble, was used as a triggering 
source for the signal processing equipment. 1 he location of each subsequent bubble was 
calculated by using the dilfercnce in time of arrival of the noise from the bubble to the 
diHercnt clcmcitts of a hydrophone array. The bubble production area can be measured 
by positioning all the bubble locations on a grid. The rate of bubble production was 
observed to behave as an exponentially decreasing function of time. 1 his was determined 
by noting the relative time of arrival of the pressuic waves from the individual oscillating 
bubbles to a hydrophone array, 'fhe bubble production density and volume rate of gas 
entrainment were calculated from the bubble population data and the measured area of 
bubble production. This technique provided information on bubble production on the 
surface, a quantity which had not been directly measured in any other previous exper- 
imental work. 



II. BACKGROUND 



Many processes of acoustical, meterological, geochemical, and biological interest are 
associated with the presence of bubbles in the near surface of the oceans. Among these 
are ambient noise, aerosol formation, surface coagulation, fractionalization of organic 
or inorganic materials, and gas exchange. The bubble populations that influence these 
phenomena are composed of bubble patches produced by breaking waves superimposed 
on a background population of bubbles. 

The background population of bubbles originates from bubble patches due to pre- 
vious breaking waves, biological activity, droplet impacts, and sediment outgassing. 
Some of these bubbles will rise and eventually burst on the surface aficeting the proc- 
esses already mentioned. Other bubbles undergo diffusion and will completely dissohe. 
.•\n expression given by Levich (1962) for the rate at which these bubbles will dissolve 
is described by the following equation: 

dn i- i- A 

-^ = -7.9SZ)3;/3u3 C 
dl 



where n is the number of moles of the gas. D is the gas diffusivity in uvj sec . u is the 
bubble velocity in m sec, a is the bubble radius in meters, and C is the Peclet number 
which desciibes diffusive mass tiansfer. 

I he vertical distribution of the bubbles beneath the surface has been studied by 
Thorpe (1986) and Farmer and Vagle (1989). Clouds of subsurface bubbles caused by 
wind waves breaking in deep water and ci eating whitecaps have been observed using a 
vcrticalh pointing narrow beam sonar and a dual beam side scanning sonar. The vertical 
distribution of bubbles is influenced by several actions: the turbulent transport from a 
source at the surfitce. the rise of the bubbles due to buoyancy, the change in bubble ra- 
dius due to the dilTusion discussed earlier, and change of hydrostatic pressure. The tur- 
bulence includes the forces from the whitecaps. the shear from the breaking wave, 
internal waves, and Langmuir circulation. Thorpe's equation representing the process is 
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Where is the rising velocity of the bubbles in m sec. N is the number of bubbles per 
unit volume, A' is the turbulent diOiision coefllcient in sec , and o is the bubble loss 
rate due to diffusion in number of bubbles per second. 

There have been many other techniques developed to measure these background 
bubbles. Johnson and Cooke (1979) and Toba (1961) utilized a photographic method 
to accomplish this task. Three images corresponding to the specular reflections of the 
light provided by three strobes served to identify bubbles and as the basis for measure- 
ment. The depth of field was established by the intersection of three dots and a high re- 
solution film allowed a 3X magnification. This method claimed identification of bubbles 
with radii as small as 17 nm. The procedure requires the water to be in a clean, clear 
chute and can identify bubble radii with only poor resolution. The abundance of debris 
and other matter would interfere with accurate identification of bubbles with a radius 
of less than 0.1 nun. Kolovayev (1976) used a bubble trap in order to concentrate the 
bubbles to allow measurements by photographic means. Breitz and Medwin (1989) de- 
veloped a floating acoustical resonator to determine the number and size of bubbles 
produced from spilling breakers at sea. 1 he Q of the resonator was changed by bubble 
presence. 'I'he system allowed measurement of bubble populations of nine radii between 
30 and 270 iim. Though accurate, this method measured bubbles at a depth of 25 cm 
below the surface and getting densities at a shallower depth may not be possible. 

With these successful results in measuring and understanding the background bub- 
bles. the ne.xt major area for study is to e.xamine the origins of the bubbles and the re- 
lationships between these bubble producing actions. 

A large proportion of bubbles in the oceans originate from the whitecaps of break- 
ing waves, llolthuijsen and llerbers (1986) conducted measurements of wave height, 
period, steepness, crest front steepness, horizontal crest asynunetry. and vertical crest 
asynunetry, in an attempt to estimate from these parameters the likelihood that a wave 
will break. They concluded that breaking waves cannot be separated from non-breaking 
waves with these factors. Thorpe and Hall (1983) arrived at several conclusions from 
observations of waves at sea: separate bubble clouds arc created by individual waves: the 
breaking crests are seldom more than several meters long; waves generally continue to 
break as they advance so the horizontal size of the bubble clouds increases with wind 
speed. Another mechanism of bubble production involves the small capillaiw waves that 
can be observed on that front face of long gravity waves even in the absence of wind. 
These capillaiw waves can break and pinch off air bubbles in the same manner as the 
breaking gravity wa\es. 
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The role of surface tension in breaking waves is also important. .Vliller (1972) per- 
formed a series of experiments on lowering the surface tension of the water with chemi- 
cals. In these experiments it was found that, as the surface tension decreased, the waves 
got higher before they broke. Gucinski (19S6) found that the persistence of visible bub- 
bles at sea is longer than that in fresh water during similar wind and sea conditions. A 
contributing factor is the presence of organic matter which surrounds the bubble inter- 
faces. Laboratory e.xperiments show an increasing coalescence with increasing ionic 
strengths. 

Since distributions of wind and wind waves can be collected by satellite sensors, the 
development of expressions for related phenomena as functions of these winds and 
waves should be useful. .Many studies have been made and great quantities of data have 
been taken of wind stress and whitecap coverage measured from a level of 10 meters. 
Toba and Koga (1968) have shown that the percentage of breaking waves ( a ). per- 
centage of whitecap coverage ( P ). and the concentration of sea salt particles ( O ), are 
related by the following term; 
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where u. is the surface friction velocity of air in m sec, v is the kinematic viscosity of 
air in sec . and o is the spectral peak frequency of the wind waves. It is now also 

possible to measure the tropospheric aerosols over the oceans with satellites. The 
method relates the upwelling visible radiants measured by the satellite to the aerosol 
ojuical thickness of the atmosphere. 1 here is a relationship between the radiance and the 
aerosol optical thickness which varies with parameters such as aerosol size, distribution, 
and refractive index. The correlation between the sources of the aerosols, such as 
breaking waves and bursting bubbles, and their concentration in the atmosphere is still 
not understood. 

The use of passive acoustics and modern signal processing allows individual bubbles 
to be examined as they arc created by the breaking wave. Medwin and Beaky (1989) 
separated the characteristics of these bubbles into four distinct categories. These 
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categories describe the bubble damping, resonance frequency, and oscillation patterns. 
With the method described in this thesis, utilizing the same equipment as Medwin and 
Beaky, the location of these individual bubbles can be deternriined from the geometry of 
the hydrophone array and the dilTerence in the time of arrival of the pressure field to the 
individual elements. From this data the exponentially decreasing temporal rate ol' bubble 
creation from an individual breaking wave can be calculated, as well as the production 
density in bubbles per unit area originating from a single breaking wave. The volume 
of gas encapsulation per unit surface area of a spilling breaker can be computed from 
the bubble density. This could be useful in studies of gas exchange at the ocean surface. 

If all the information and methods described can be integrated and successfully uti- 
lized it will be possible to obtain a complete perspective on the generation of bubbles 
at sea. The background bubbles are understood and can be measured; the main source 
of bubble creation, particularly whitecap coverage from breaking waves, can be esti- 
mated; the physical and chemical makeup of the water can be examined; and now with 
the use of the method described in this thesis it is possible to predict the rate of bubble 
production from the whitecaps. This information, when combined, will further aid in our 
understanding the efl'ect of air sea interaction on our environment. 
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111. FACILITIES AND APPARATUS 



The lahoraton" facilities include a 57x4x4 foot wave tunnel emptying into a 
10x10x10 foot anechoic tank. A motor-driven reciprocating wedge generates the waves 
at the beginning of the tunnel. The waves break approximately three feet after entering 
the tank. The exact position of where the waves will break was not controlled. The tank 
sides and bottoms arc lined with Redwood pilings. .After spilling, the waves are absorbed 
by a "beach" consisting of an aluminum shavings wedge which minimizes the wave re- 
flection at the tank boundarv'. 

Medwin and Beaky (1989) conducted tests in the tank which demonstrated at least 
a 15 decibel dilTerencc between a signal and the strongest reflection from the tank walls 
at frequencies under 5 kHz. and at least 20 dB at frequencies above 5 kHz. when the 
hydrophone is 24 cm from the bubble source. The signal to reverberant noise ratio is 
better for shallower hydrophones. 

Two vertical arrays of 2 onmidiicctional hydrophones with the shallow hydrophone 
at a depth of 15 cm and the deep ht drophone at a depth of 31 cm are situated in the 
approximate area where the majority of the waves break as depicted in figure (1). The 
output signals arc input into an Ithaco 1250 voltage amplifier bandpass filter with the 
bandpass frequencies set between 300 Ilz and 100 kHz. 'I he amplifier output splits into 
a Hewlett-Packard 3561A signal analyzer and also into an IB.M XT personal computer 
with an R. C. Electronics "comptitcrscope" analog to digital converter, fhe signal ana- 
lyzer and the "computerscope" allow the operator to select triggering level, 'fhe signal 
analyzer triggers at an input level of 0.837 Volts which, for an amplification of 5000. 
corresponds to an acoustic pressure of 1.4 Pascals. 

I his "computerscope" analog to digial conversion software package samples the 
data with a sampling period of 1 microsecond with one input. 2 microseconds with two 
inputs, and accordingly. 4 microseconds with four inputs with a 12 bit amplitude resol- 
ution. The sampled signal is displayed both graphically and nunicrically with cursor 
movement. The buffer is 64 milliseconds long with the minimum sampling period of I 
microsecond and can be lengthened by a corresponding increase in the sampling period. 
Sections of the buffer can be "expanded" or "compressed" in time and the amplitude of 
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the signal display can also be adjusted. This enables the user to manipulate the graphical 
presentation to obtain a signal that can be accurately analyzed and also to adjust the 
sampling rate to maximize the signal resolution. 
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Hydrophone 




Figure 1. Equipment setup. 
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IV. EXPERIMENT 



A. BREAKER SPECTRUM 

During World War II many observations of the ambient noise of the oceans showed 
that "wind noise" is the primar\’ source of the noise in the frequency range between 500 
llz and 20 kllz. The famous Knudsen spectra broke down this large collection of data 
in a manner intended to demonstrate the wind speed dependence of the noise level. This 
spectrum is reproduced from Urick (19S3) as figure (2). The possible processes, as given 
by Urick, which produced this noise were hypothesized to be "crash noises of breaking 
waves", "flow noise" produced by wind moving over the water, "cavitation noise", and 
"wave generating actions" of the wind. The paper by Medwin and Beaky (19S9) proves 
that the catastrophically created, damped bubbles, by themselves, explain the Knudsen 
spectra. The slope of the noise spectrum in this frequency band is approximately -5 to 
-6 dB per octave or -17 dB per decade with the level increasing with wind speed. The 
noise in the frequencies outside this band was believed by Urick to originate from 
sources not directly dependant on wind speed such as shipping below 500 llz and ther- 
mal noise above 50 kl Iz. 

In order to determine the noise spectrum of the waves breaking in the tank it was 
first necessary to ensure that the wave generating equipment was not acting as an addi- 
tional noise source. Figure (3) shows the background ambient noise of the tank with no 
equipment in operation plotted with a dashed line. The ambient noise in the tank with 
the plunger in operation but no waves breaking is plotted with a solid line. Both plots 
were generated on the IIP 3561,\ spectrum analyzer using a time average of 3UU input 
signals, each of 5U0 milliseconds duration. The signal level was computed for Irequencies 
between 500 Mz and 50 kl lz. 

The breaker spectrum was also obtained using the IIP 3561 A with the input from 
an onmidirectional hydrophone at a depth of 24 cm. The signal was divided into 40 se- 
quential records, each of 20 millisecond duration. Each 20 ntillisccond band was Fourier 
transformed to provide the power for each 50 llz block of the bulfer.'*^ fhe energy in 
these bands was summed only if the energy was 3 dB above the noted ambient noise level 
of the tank with no activity. This was done to ensure that the energy originated only 
from breaking waves in the tank. The breaking of a wave can usually be identified by 
a relatively high amplitude low frequency bubble which served as the triggering source 
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for the signal analyzer. The frequencies that were observed were restricted by the spec- 
trum analyzer to between 400 llz and 20 kHz. A total of six breaking waves was meas- 
ured and the levels averaged. The results are displayed in figure (4). The slope of the 
spectrum is similar to the Knudsen spectra at approximately -17 dB per decade. The 
noise level at a frequency of 1 kHz is approximately 82 dB re liuPa-jlJz which is some- 
what greater than that of the Knudsen spectra at that frequency at high wind speeds. 
This is due to the laboratory measurements being made only when and where bubbles 
are being produced; a comparable situation at sea would require breaking waves to exist 
ever>"where at all times as discussed by .Medwin and Beaky (1989). The low noise level 
of figure (3) confirmed that the only significant sources that contribute to the sound 
levels above the ambient room noise are breakers produced by wave actions and not 
mechanical equipment. It should be observed that this phenomenon was observed tinder 
the condition of no wind over the tank. This confirms the suspicion that the actual 
source ol' the ambient noise are the breaking waves and not the wind. It was the fasci- 
nating fact that the spectrum in the tank is \ery similar to the expected spectrum in the 
oceans that convinced the author that the data produced by this study may be of future 
use to actual operations performed at sea. The method of using wind speed as a reference 
is useful, however, because the wind speed can be more dependably measured than the 
sea state. An interesting thing to note about the breaker spectrum is that below 1 kHz 
the level remains constant to the minimum frequency measured at 4U0 Hz. 

.M. S. Longuet-1 liggins picdicted that certain frequencies would be predominant 
from bubbles produced by breaking waves. 1 Iiese frequencies originate from modal dis- 
tortions in the shape of the bubble during its creation. '1 he frequencies were stated in a 
submitted paper. Longuet-1 liggins (1989). to be at 1.38, \.09. 2.10. 2.60, 3.32, 4.31, .3.70, 
7.94. 1 1.4, and 16.8 kHz. Ihesc frequencies arc designated by triangles on figure ^4). 'Hie 
laboratory spectrum sliowed a peak only at two ol' these frequencies, 1.69 and 2.6>0 kHz. 
Several of the other frequencies correspond to the position of a minimum on the spec- 
trum. T his docs not allow a determination of the significance of this phenomenon. 
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Figure 2. Average deepwater ambient noise spectra. (Urick) 
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Figure 3. Ancchoic tank background noise spectra: Solid line, plunger operating 
but no breakers. Dashed line, plunger off. 
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Figure 4. Average noise spectrum from six breaking waves: Triangles mark the 
frequency peaks predicted by Longuet-Higgins. (1989) 
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B. BUBBLE IDENTIFICATION 



It is; possible to positively identify individual bubbles within the sound originating 
from the breaking of a wave. The equipment utilized is a hydrophone, a preamplifier, 
and an IBM XT personal computer with an R.C. Electronics Computcrscope analog to 
digital converter. The converter samples the signal and allows detailed viewing of the 
desired segments of the captured noise. 

^^'hen a wave breaks or "spills" air is entrained, bubbles are formed, and the bubbles 
oscillate to become sources of noise. In eveiw wave observed with the signal processing 
equipment, the moment when the wave breaks can be identified by a large amplitude, 
low frequency bubble, with sulEicient energy to be used as the triggering source for the 
signal analyzing equipment. 

This equipment has already been used by .Medwin and Beaky (19S9) to characterize 
the bubbles from breaking waves by identifying the damping and modulation aspects 
most conunonly observed. Four distinct bubble scenarios were classilied: spherical bub- 
bles with either a constant decay rate or with two dilTcrcnt decay rates, a damped oscil- 
lation bubble appearing to shed a much smaller bubble, a near surface bubble which has 
an increasing amplitude as it moves away from the surface due to an increase in the 
dipole axis, and a bubble undergoing amplitude modulation possibly due to another 
bubble in close proximity with a slightly diflcrent frequency. 

The resonance frequency of the bubbles was converted to bubble radius using the 
Ibllowing relationship from Clay and Medwin (1977); 



where f is the resonance frequency of the bubble in Ilz, a is the bubble radius in meters, 
y is the ratio of specific heats of the gas in the bubble, P is the hydrostatic pressure in 
.V/'/nk and p is the density of the water in kglnt^. 

’1 he large amplitude initial bubble created by a breaking wave can be seen in figure 
(5). This is the first 120 milliseconds of a 5o0 millisecond bulTer. In this printout the 
upper and lower signals are from a hydrophone at 16 cm depth and a hydrophone at 31 
cm depth respectively. The output signals arc ISO degrees out of phase because of the 
amplifier output phase dilference. Ihe confirmation that it is indeed a bubble can be 
made from the damping relationship from Medwin and Beaky (1989) 
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where t is the decay time in seconds to 1 e of its amplitude, fis the resonance frequency 
in Hertz, and S is the damping constant. The measured frequency of the bubble in figure 
(5) is 505 Hz and the measured l, e decay time is 30.9 milliseconds. From the table on 
page 199 of Clay and Medwin (1977) the damping constant for this frequency should 
be approximately 0.022. If placed in the above equation this damping constant yields a 
I e decay time of 28.6 milliseconds-, a difference from the measured value of 7 percent. 
This is close enough to the theoretical value to ascertain that the bubbles are behaving 
as expected. If the presentation is expanded and viewed with more detail as in figure (6), 
additional later signals from higher frequency, smaller amplitude bubbles can be seen 
superimposed on the low frequency bubble signal. Finally, the presentation can be ex- 
panded even further, figure (7), to enable the accurate identification of a smaller bubble 
with frequency 16.6 kl Iz (196 micron radius) and accurate measurement of the difference 
in time of arrival of the signal to the individual hydrophones. 

It can be seen in the above figures that the individual bubbles can be identified with 
enough resolution to obtain accurate information about bubble frequency, amplitude, 
and relative time of arrival at the hydrophones. 1 he signal data was recorded only if a 
bubble could be confidently identified above any noise which may have been received. 
^\'ith this equipment accurate identification could be ensured if the signal level was at 
least 3 decibels above the measured ambient noise level of 37 dB re l/.iPa^lIIz at the 
hydrophones. The resulting bubble count is thereby a conservative measure of the 
number of bubbles and would not include any bubbles with received signal less than the 
limit of 37 dB. 
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6 uS 24.6 nS/Div 126,066 uS 



Figure 5. Large triggering bubble from breaking wave; Predominantly large 
amplitude and low frequency bubble (505 Hz). Signal presentation is 120 
ms. Outputs from two hydrophones shown. One output is invened by the 
amplifier. Sampling period is 15 microseconds. 
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Figure 6. Additional smaller bubble identification: Several bubbles with higher 
frequency can be seen superimposed on the initial bubble signal. This is 
15 ms of the 120 ms plotted in figure 5. Sampling period is 15 microsec- 
onds. 
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Figure 7. Small bubble identincation: Frequency is 16.6 kHz. This is a 1.875 
ms section of the signal shown in figure 5. Sampling period is 15 micro- 
seconds. 
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C BUBBLE LOCATION 



The location of the bubble can be determined using the geometry of the hydrophone 
positions and the relationships developed from the difference in time of arrival of the 
sound pressure wave to the hydrophones. In this procedure it was assumed that the 
bubbles were created and remain approximately at the surface of the water. The argu- 
ment for this assumption follows: the bubbles were produced from the action of the 
breaking wave entraining pockets of air. This catastrophic force was the driving mech- 
anism causing the oscillation. It was observed that the bubble oscillations only lasted for 
a short duration, of the order of tens of milliseconds, after the bubble formation. In this 
short period of action the bubble could not physically travel any significant distance. 

A second assumption is that the surface of the water remained horizontal, ignoring 
the wave height. The error due to this assumption will be examined shortly. 

\Mien two hydrophones are placed in a vertical array the sound pressure wave from 
•a bubble at the surface will arrive at the hydrophones at different times. This difl'erence 
in time of arrival will be a function of the distance the wave must travel through the 
water. 

Bubbles produced near the surface of the water act as dipole sources. The dipole 
consists of the bubble and an imaginary source the same distance above the surface as 
the bubble is below, with a distance L separating this image and the bubble. With the 
assumption that kL < 1. the dipole radiation pattern will be defined by the following 
equation from Morse and Ingard (196S): 

P = ( )pcD cos(f?)( 1 + 

where k is the wave number in radians meter. R is the range in meters Ifom the bubble 
to the hydrophone, pc is the acoustic impedance of water in Pa sec m. D = VL where 
f^=4,-TU-L' is the volume per time. U is the radial velocity amplitude. 0 is the angle be- 
tween the axis of the dipole and the line adjoining the bubble and the hydrophone, and 
w is the angular frequency of the oscillating bubble in radians second. The magnitude 
of the pressure at hydrophone 1 can be described as 

\}\\ = ik^pcD cos( )][ 1 + ]T. 



Using the relationship 
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a + ih = {a^ + b^) 2 



so that the term 



{ 1 + — ^ ^ = n 4- — ^ — ) 2 
kR ’ ^ /^. 2;^2 ) 



and 



0 = tan ’(^) = tan '(“tV). 



The real part of the pressure at hydrophone 1 can now be described by the equation; 






or restated 



A, = I A, I cos(u7, - kR^ + tan '( -pj— )). 

f\ /< j 

Identically, the real part of the pressure at hydrophone 2 can be written: 

P, = I Aj 1 cos(r.7, - A- AN + tair'( )). 



Comparing the phases for the times that the peaks occur, 

u/j — AA, + tan~'( ~ r \ - ) = n~ 

k /\i 



and 



u72 — a Ai + tan ’( ! ) = tire 

k lx-) 



and subtracting yields: 

u -(/2 - - A(A 2 - A,) + tan-'( ^ ) - tan-'( -^ ) = 0. 

T he equation describing the diflerence in time of arrival can be written as; 
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^,=(- 



R. - R, 



) + (■ 



• w ' 

T ) tan ( 



Inf 



kR, 



^ ^ 2nf ^ ^ ^ 



where /?, and R 2 are the distances from the bubble to the hydrophones in meters, k is the 
wave number in radians meter, c is the speed of sound in meters second, and f is the 
resonance frequency of the bubble in Hertz. 

The geometry of the hydrophones is depicted in figure (8). The following relation- 
ships should be noted: 

/?, = (/if -y .x^)2 
= (/'2 + 

where and lu are the depths of the hydrophones in meters and x is the horizontal dis- 
tance on the surface of the water above the hydrophone array in meters. If the position 
on the surface above the h} drophones is treated as an origin, x is in fact the radius of a 
circle. If these equations for /^, and R 2 are placed in the previous equation for time of 
arrival difference, this radius, x, will remain the only unknown. Using an iterative pro- 
gram. the computer will easily determine the solution. 

The position of the bubble on the surface can be determined by using two vertical 
arrays of hydrophones. The bubble w ill be located somewhere on the circles surrounding 
the arrays. These circles will intersect at two points, in a symmetrical manner about the 
horizontal line separating the arrays as in figure (9). The ambiguity between the two 
intersections is unimportant to the objectives of the experiment. 

The fact that the surface of the water is not constantly horizontal, but in reality 
consists of various wave heights, causes a change in the time of arrival dilference be- 
tween the hydropliOnes of an array. If the bubble is located on a war e crest or trough 
the distance the pressuie wave will ha\ e to travel can be described by: 

/?, =(f/i, -b A/i)^ + .v^)i 

R2 = {{I12 + Mif -t- ■*' 



where A/i is the vertical difTerence in bubble position from the horizontal surlace of the 
water, 'fhis difference in time of arrival caused by the wave amplitude is shown in figure 
(U>). It will be noted that the maximum amplitude of the waves in the tank is 



21 



approximately 3 cm. At this wave height the error in time of arrival is only 5 microsec- 
onds. Because of the fact that the wave height follows a Gaussian distribution, the wave 
height at the bubble location will likely be on the order of 1 cm and the error in time 
would be 2 microseconds with a corresponding error in bubble location of an insignil- 
icant 2 or 3 millimeters. 

Two vertical arrays, each including a shallow hydrophone at 16 cm depth, and a 
deep hydrophone at 31 cm depth, were placed in the tank in the area where most waves 
break. The perpendicular line between the hydrophone arrays was set at approximately 
a 30 degree angle to the direction of the wave propogation in order to ensure geometric 
efficiency in the calculation of bubble positions. The buHcr length was set to 100 milli- 
seconds with a corresponding sampling period of 6 microseconds (spatial resolution, 9 
nun). This small sampling period was considered necessary in order to measure the time 
of arrival dilTerence as accurately as possible. The output went through the preamplifiers 
set at a gain of 500o and acted as the trigger to the computerscope analog to digital 
converter. 

'1 he dilference in time of arrival to the two vertical pairs of hydrophones, the reso- 
nance frequency of the bubble, and the time that the bubble began oscillating relative 
to the breaking of the wave, were recorded in the first 100 milliseconds for each bubble 
created by a breaking wave. This information was used with the iterative computer 
program to plot the locations of the individual bubbles. Tigure (11) is the result of one 
breaking wave showing the synunetrical distribution of the bubbles created in the first 
lot) milliseconds after breaking. Ihe synunetrical locations are caused by the two inter- 
sections of the bubble radii circles about the axes of the hydrophones. .As long as the 
wave does not break close to the hydrophones it is not necessary to plot both sets of 
location positions because it is the relative positions of the bubbles that is necessary for 
density calculations. I'igures (12) and (13) show the bubble locations with only the pos- 
itive positions recorded. 

1 he position plots demonstrate the bubble patches or "hot spots" olten noted in 
acoustical measurements at sea. The plots were divided by a grid consisting of squares 
5 cm on a side. These blocks of 25 square centimeters were used to calculate the area 
over which bubbles were produced. This large area assigned to any bubble counters and 
compensates for any error in the calculation of bubble position. l:ach square with a 
bubble in residence was counted. 1 he average total area of bubble production was 345 
square centimeters per breaking wave with a standard deviation of 42 square centime- 
ters. It was assumed that all bubbles created by a breaker originate in such a calculated 



area. It was not possible to expand the time of observation in order to confirm this as- 
sumption for later bubbles because the sampling period would also have become longer 
and the resolution of the time of arrival dilference would not have been sufficient for 
accurate bubble location. Another method was to form a rectangle around the bubbles 
and compute the production area. The average production area calculated with this 
method was 320 square centimeters with a standard deviation of 96 square centimeters. 
The rectangle is marked by a solid line around the bubbles in figures (12) and (13). The 
two methods compare within 10 percent for the average production areas. 

Another result of the bubble position plots was the ability to label the order in which 
the bubbles were created by the breaking wave. Figures (12) and (13) also show these 
relative positions of the bubbles from two waves each marked with the order in which 
they were produced. The order of the positions does not demonstrate that the location 
of the bubble depends on the relative time of its creation with respect to the breaking 
of the wave, 'fhe first bubble created was not necessarily the closest to the position 
where the wave breaks. Also the last bubble created was not necessarily located the far- 
thest along the direction of wave propagation from the position where the wave began 
breaking. This fact suggests that the majority of the bubbles created were located within 
the area calculated even though this area was calculated from bubbles originating in the 
first lUU milliseconds of breaking. 

Bubbles observed within the bubble location plots were ocassionally positioned in 
roughly linear patterns not necessarily perpendicular to the direction of wave 
propogation. The bubbles within these patterns all originate consecutively in short 
widows of time. It is believed that these bubbles are created by and along the small 
capillary waves which sometimes are seen superimposed on the face of a main wave. This 
evidence of capillary waves is examined in Appendix A. 
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Figure 8. Bubble location geometry. 



24 



Possible 

Bubble 

Positions 




Figure 9. Bubble location technique: Intersection of the two circles of radii x ^ 
and X , around the hydrophone axes marks the possible location of the 
bubbles. 
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Hrror in time of arrival difference (j(s) 




Horizontal surface distance from bubble to array axis (m) 



Figure 10. Timing error due to wave height: Error in time of arrival caused by 

the additional distance traveled by the noise of a bubble on the surface 
of the water due to local wave height above mean surface. 
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Figure 11. Distribution of bubble locations; Distribution is symmetric about the 
line joining the two vertical arrays. Bubbles are from the first 100 ms 
after the breaking of the wave. 
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Figure 12. Bubble location plot: Bubble locations marked with the order of bubble 
creation. Same breaker as figure 11. Bubbles are from the first 100 ms 
after the breaking of the wave. 
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Figure 13. Bubble location plot: Bubble locations marked with the order of bubble 
creation. Bubbles are from the first 100 ms after the breaking of the wave. 
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D. RATE OF BUBBLE PRODUCTION 



It was observed from the use of the Hewlett-Packard 3561A Spectrum Analyzer that 
noise is produced from a breaking wave for a period of 2 to 3 seconds after the first 
sound of a breaking wave is heard. The digital processing equipment used to identify, 
locate, and count bubbles lias a relatively small memory (64k bites). This does not per- 
mit the sound to be captured for a long enough period of time to record all bubbles 
produced and simultaneously have a sampling period small enogh to be of use in the 
identification of bubble resonance frequencies. Therefore it was thus necessaiy to deter- 
mine the rate at which bubbles are produced over time in order to determine an e.\- 
pression defining the bubble production rate relative to the time at which the wave 
breaks. The measurements could then be made for shorter periods of time and then e.K- 
trapolatcd to determine the total number of bubbles produced. 

Two omnidirectional hydrophones were placed in a vertical array in the tank at 
depths of 12 and 24 centimeters. The buirer length for two hydrophones was expanded 
to 500 milliseconds by using a sampling period of 15 microseconds. Through exper- 
imentation and consideration of the Nyquist relationship it was determined that this was 
the maximum sampling period that would allow dependable identification of bubbles 
with resonance frequencies less than 33 kHz. The trigger level was adjusted to a level 
sulficiently high to ensure the sampling equipment would begin capture only when a 
wave was breaking almost directly above the hydrophones. 1 his was done to guarantee 
the obsen ation of bubbles with a small amplitude which may not have been seen from 
a longer distance and at a greater angle from the dipole. The short range from the 
hydrophones can be confirmed using the method previously described. The input signal 
passed through the preamplifier with the gain set on 5000. A total of 10 breaking waves 
were observed. 

The average number of bubbles of all radii produced per breaker, within time incre- 
ments of 10 milliseconds, is plotted against the time after the breaking of the wave in 
figure (14). Additionally, the production rates for bubbles with bubble radii broken down 
into bands of greater than 2.2 mm (resonance frequency less than 150u Hz), between 2.2 
mm ( 1500 Hz) and U. 10 mm (20 kl Iz). and less than 0.16 mm (greater than 20 kHz), are 
depicted in figures (15). (16). and (17) respectively. These bands are chosen to accentuate 
the bubble production behavior in the Knudsen sea noise spectrum fretjucncy regime 
generally taken to be between 500 Hz and 20 kHz. 

T'he production spectrum of the two radii bands. 0.16 to 2.2 mm and smaller than 
0.16 mm as seen in figures (16) and (H). suggests that bubbles in these bands are 
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produced witli an exponential rate of decrease. The equations which approximate tlie 
decay rates for these two bands arc: 

r/,V A -! 1 1 
— — = 0.31 le 
di 



and 



di 

respectively. This is probably also true for the largest bubbles but not enough data is 
available to show it. 

The curve which best approximates the exponentially declining total production rate 
of figure (14) is: 



r/.V -> < -> -.0070/ 

dl 

where dN dt is the number of bubbles produced per millisecond and t is time after the 
breaking of the wave in milliseconds. The expression stated above can be integrated to 
describe the bubble production rate out to specified times. \Mien the integration was 
performed to infinity it was found that less than 0.1 percent of the bubbles would origi- 
nate after the first second. About 97 percent of the bubbles were formed in the first 500 
milliseconds after the wave breaks. 
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Average number of bubbles produced 




Time after breaking of wave (milliseconds) 



Figure 14. Average number of bubbles produced per laboratory breaker; (all 
sizes). The bubbles are counted within 10 ms increments. 
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Figure 15. 




Time after breaking of wave (milliseconds) 



Average number of large bubbles per laboratory breaker: 

*C 

2.2 mm). The bubbles are counted within 10 ms increments. 



(radius > 
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Average number of bubbles produced 




Time after breaking of wave (milliseconds) 



Figure 16. Average number of mid-size bubbles per laboratory breaker: (0.16 
mm< radius< 2.2 mm). The bubbles are counted within 10 ms increments. 
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Figure 17. Average number of small bubbles per laborat 

0.16 mm). The bubbles are counted within 10 n 



V breaker: (radius < 
.ncrements. 
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E. BUBBLE PRODUCTION DENSITY 



The rale at which bubbles will be produced is related to the foam coverage. With 
this knowledge researchers will be better able to understand the phenomena which affect 
acoustics, weather, aerosol production, and any other sea surface interactions. The two 
major requirements for densil\' calculations have already been covered: the bubble cre- 
ation rate relative to the start of the breaking of the wave, and the area over which these 
bubbles are being produced. 

All bubbles in figures (15) - (17) with equal radii were summed and the average 
number of bubbles per breaker with that radius was computed by dividing by the num- 
ber of waves obser\ ed. In an earlier experiment it was found that only 3 percent more 
bubbles originated in the second 500 milliseconds after the breaking of the wave. Because 
the measurements were again taken over a period of 500 milliseconds the number of 
bubbles observed was thus increased by 3 percent to give the number of bubbles in the 
first second. This number is now defined as the "total" number of bubbles. In an addi- 
tional earlier experiment it was found that the average area of bubble production 
spanned 3-15 square centimeters or equis'alently 0.0345 square meters. The "total" 
number of bubbles was divided by 0.0345 to give the resulting average bubble production 
density in bubbles per square meter for each bubble radius observed. Thus the bubble 
production will be expressed in bubbles per square meter and assumes negligible num- 
bers of bubbles are Ibnned beyond the first second after the wave breaks. 

Bubbles wore observed with finite numbers of different size ladii. If an infinite 
number of waves were measured, eventually bubbles would be observed having every 
possible radius up to the maximum gas volume that can physically be encapsulated. To 
account for this limiiation in the data recorded, the span in bubble radius in microns 
was calculated to include half the bubble radius size between each of the subsequent 
smaller and larger bubble observations. The number of bubbles per square meter was 
divided by this radius span to give the bubble density in bubbles per square meter per 
micron of bubble radius inciement. 

The bubble production density is plotted on semilog paper with the radius on a 
logarithmic scale in figure (18). and both bubble radius and density on a logarithmic 
scale in figure (19). It can be seen that the most commonly produced bubble was 0.15 
mm in radius produced at the approximate density of 5.9 bubbles per square meter per 
micron radius increment. 1 he smallest bubble radius measured was 0.04S mm which has 
a corresponding resonance frequency of 66 kllz. The largest bubble was of radius 7.4 
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mm with a corresponding resonance frequency of -140 Ilz. It will be assumed that still 
larger bubbles would not have been found even if more breakers had been observed. 

Once the production density of bubbles for each radius has been calculated the vol- 
ume of gas encapsulated from each known bubble radius can be computed. Figure (20) 
is a sumnaary of cubic centimeters of gas volume encapsulated per square meter of 
breaking wave for bubbles of various radii. Although it was seen in figures (18) and (19) 
that the majority of bubbles produced were less than 0.5 imn in radius, the contribution 
to total gas entrainment from these small bubbles is nhnimal. Assunting that there are 
no bubbles larger than the largest measured in 10 breakers which had a radius of 7.4 
mm. the individual contributions from each bubble radius can be summed to give a 
volume rate of gas encapsulation which averages 2 } cubic centimeters of gas 
encapsulated jicr square meter of breaking wave in the laboratorv'. More than half of 
this total volume came from bubbles with radius greater than 3 millimeters. 
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Figure 18. Average surface density of bubbles produced. 
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Figure 19. Average surface density of bubbles produced. 
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VOLUME (cunic CM/Sq.M) 
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Figure 20. Volume of air encapsulation per surface area. 
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V. COMPARISON WITH OTHER DETERMINATIONS 



Wc know of only one other attempt to measure tlie bubble production in a spilling 
breaker, that was by Toba (1961). Contrary to our breakers being produced by an os- 
cillating wedge with no wind. Toba used wind of difTerent speeds and fetch in a Hume. 
Toba's method considered that an equilibrium existed between the bubbles entering from 
periodic breaking at the water surfitce and those rising to the surface. Ilis photographs 
were made of the volume density of bubbles, not the surface production density of the 
bubbles which is our goal. 1 le used the time between breakers to infer the bubble 
encapsulation rate using the relationship 



B = {-^) — 

^ ^ 100 ’ T 

where /? is the bubble encapsulation rate in bubbles square cm second. is the per- 
centage of crests actually entraining bubbles. 0 is the total number of bubbles of a par- 
ticular size in a water column of unit base area, and T is the period of the main surface 
wave. 1 he shape of Toba's encapsulation spectrum, figures (21) and (22). was very sim- 
ilar to that in this experiment with a sharp peak at a single bubble radius. However the 
peak of Toba's bubble production was at approximately 0.38 mm in radius at a rate of 
approximately 0.2 bubbles per square centimeter per second for a wind speed 12.1 m sec. 
Ilis photographic measurements included 6 different sizes of bubbles with radii from 0.15 
mm to 6.0 mm. Our method observed 53 bubble radii from 0.048 mm to 7.4 mm. His 
smallest radius bubble was our most common bubble. The production rate inferred by 
Toba for 12.1 m sec wind was of the same order as the rate measured directly acous- 
tically. 

'I here have been extensive measurements of the foam coverage at sea using aerial 
photography and satellite imaging. The foam coverage is a measurement of residual 
bubbles including those convccted upward from below and lacking those that have been 
collected from the surface to lower deinhs or that have burst at the surface, 'fhus the 
area of coverage can not be used directly with the densities measured in this experiment 
to calculate the bubble populations beneath the ocean surface. Further experimentation 
is required in order to understand exactly what is being measured from the foam of 
breaking waves. 
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Figure 21. Rate of formation of air bubbles in wind waves (Toba). 
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Figure 22. Rate of formation of air bubbles in wind waves (Toba): plotted per 
equal log (bubble size) interval. 
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VI. CONCLUSIONS 



The method described throughout this thesis achieved its goal of obtaining the 
bubble production density from plunger generated spilling breakers. The technique used 
was to passively observe the amplitude and frequency of the randomly pulsed sound 
emitted by the breakers. The results give more details of the characteristics of bubble 
production than previous experimental techniques used for obtaining surface bubble 
densities. 

Bubbles were positively identified throughout a radius range of 0.04S to 7.4 milli- 
meters. Most of the bubbles originate in the first .500 milliseconds after the breaking of 
the wave though sound from bubbles has been recorded for a period of 2 to 3 seconds 
after breaking. Ihe bubbles were produced with an exponentially decreasing rate fol- 
lowing a decay constant of-0.007t with time measured in milliseconds after the breaking 
of the wave. The average area of bubble production from breaking waves in the tank 
was 345 square centimeters. 1 he most commonly produced bubble had a radius of 0.15 
millimeters and was entrained at a density of 5.9 bubbles per square meter per micron 
of bubble radius increment. This bubble radius is smaller than the minimum bubble ra- 
dius measured in other experiments using photographic techniques to obtain bubble 
volume density. The greater resolution using acoustics instead of photography has al- 
lowed 10 times the resolution in bubble radii identification for bubbles with radii be- 
tween 0.25 and 2.0 millimeters. The average encapsulation of gas was 23 cubic 
centimeters per square meter of breaking wave. 

Using the Knudsen sound radiated from the breaking waves as our sole criterion, 
whitecaps from spilling ocean waves have little physical dilTerencc from the breaking 
waves produced in the laboratory. If this is the case, the measurements reported here 
are useable for the ocean. The use of this passive acoustical technique in actual sea 
measurements may also be possible. If these laboratory data are confirmed by actual 
measurements of surface bubble production density at sea, and after further studies of 
bubble dissipation in a turbulent ocean, satellite measurements of the foam coverage 
from whitecaps could yield an estimate of the rate of bubble production in the ocean. 
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APPENDIX A. EVIDENCE OF CAPILLARY WAVES 

Upon examination of the bubble location plots it was discovered that in some of the 
plots the bubbles appeared to be located along lines perpedicular to the wave propa- 
gation direction (fig 12) and in other cases the lines are not perpendicular to the direction 
of wave propagation. These lines appeared with varying angles ranging from 30 to 90 
degrees away from the propagation direction. It was hypothesized that these may have 
been capillary waves forming on the face of the main wave. A purely arbitrary' reference 
time of 10 milliseconds was chosen as the time within which all of the bubbles on the line 
must have originated. Tigtire (23) is an example of a capillan.' line identification. The 
bubbles were captured with a sampling period of 6 microseconds during the first 100 
milliseconds after the wave broke. The order of bubble formation is marked on the fig- 
ure and the following are the times of the beginning of the individual bubble oscillations 
after the breakinc of the wave: 



: number 


Origin time (ms) 


1 


2.92 


2 


5.88 


J 


10.56 


4 


11.22 


5 


12.44 


6 


18.74 


7 


20.54 


S 


24.99 


9 


31.15 


10 


33.66 


11 


42.37 


12 


42.49 


13 


49.03 


14 


56.52 


15 


69.30 


16 


75.23 


17 


78.65 
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IS 

19 

20 



80.45 

81.83 

93.94 



The capillan' lines on figure (23) are identified by bubbles 6 through 8 and by bubbles 
17 through 19. Figures (24) and (25) are photographs of breaking waves in which the 
capillary waves may be clearly observed. 
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Figure 23. Possible capillary wave positions b\ bubble location: Lines are from 
bubbles 6-8 and 17-19. 
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Figure 24. Photograph of capillary waves. 
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Figure 25. Photograph of capillary waves. 
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APPENDIX B. VOLTAGE AMPLIFIERS 

When receiving signals from many hydrophones simultaneously, it was necessaiy to 
develop a substitute for the limited number of Ithacc 1201 Preamplifiers. A series of 
preamplifiers was built using two 356 operational amplifiers in series with the necessar>' 
feedback elements. Figure (26) is the diagram of the amplifier circuit. Figure (27) shows 
the operational characteristics of the amplifier. 
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Figure 26. Amplifier circuit diagram. 
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Figure 27. Operational characteristics of the amplifier. 
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